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Abstract 

This paper examines convergence and stability of the two classes of theta-Milstein schemes 
for stochastic differential equations (SDEs) with non-global Lipschitz continuous coefficients: 
the split-step theta-Milstein (SSTM) scheme and the stochastic theta-Milstein (STM) scheme. 
For 9 G [1/2,1], this paper concludes that the two classes of theta-Milstein schemes converge 
strongly to the exact solution with the order 1. For 6 € [0,1/2], under the additional linear 
growth condition for the drift coefficient, these two classes of the theta-Milstein schemes are also 
strongly convergent with the standard order. This paper also investigates exponential mean- 
square stability of these two classes of the theta-Milstein schemes. For 9 G (1/2,1], these two 
theta-Milstein schemes can share the exponential mean-square stability of the exact solution. 
For 9 G [0,1/2], similar to the convergence, under the additional linear growth condition, these 
two theta-Milstein schemes can also reproduce the exponential mean-square stability of the 
exact solution. 
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1 Introduction 


During the last decades, SDEs have become increasingly important tools to describe the real world 
since stochastic models have wide applications in biological systems, neural network, financial 
engineering and wireless communications. Most SDEs cannot be solved explicitly, so numerical 
approximations become important tools to study stochastic models. When a numerical scheme is 
put forward, it is crucial that this numerical scheme can converge to the exact solution. Moreover, 
it also need to describe the asymptotic properties of the exact solution such as boundedness and 
stability. 

Most of the existing convergence theory for numerical methods of SDEs requires the global 
Lipschitz condition (see m)- However, many well-known stochastic systems do not satisfy the 
global Lipschitz condition, for example, the stochastic Duffing-van der Pol oscillator m, stochastic 
Lorenz equation EE], experimental psychology model [5], stochastic Ginzburg-Landau equation 
[HE], stochastic Lotka-Volterra equations mum and volatility processes mu and so on. For 
some SDEs without the global Lipschitz condition, the classical explicit numerical schemes may 
not converge to the exact solution in the strong mean-square sense (for example Euler-Maruyama 
(EM) scheme and Milstein scheme, please see [Hl[9|). Hence, numerical approximations for SDEs 
without the global Lipschitz condition have received more and more attention in recent years. 

For the SDEs with the one-sided Lipschitz condition (which is weaker than the global Lipschitz 
condition) on the drift term and the global Lipschitz condition on the diffusion term, Hu |10] 
examined convergence of the backward Euler-Maruyama (BEM) scheme and obtained that the 
optimal rate of convergence is 0.5. Under an additional polynomial condition, Higham et al. unna 
and Bastani et al. m proved that BEM and split-step BEM schemes converge strongly to the exact 
solution with the optimal rate 0.5. Mao and Szpruch HD showed that under a dissipative condition 
on the drift coefficient and the super linear growth condition on the diffusion coefficient, the BEM 
scheme is convergent with strong order of a half. Under a monotone condition, strong convergence 
of BEM scheme and theta-EM scheme were investigated by [151116] . Recently, there are various 
explicit Euler schemes with one half order were proposed to approximate the SDEs with non-global 
Lipschitz coefficients, see [IIH2ii| and the references therein. However, it is still an important and 
difficult work to look for more appropriate conditions to obtain the strong convergence rate of the 
numerical schemes for SDEs without the global Lipschitz condition. 

The explicit Milstein scheme for SDEs developed by Milstein [2] can obtain a strong order of 
convergence higher than Euler-type schemes. Hu et al. [21] extended this scheme to solve SDEs with 
time delay. Implicit Milstein schemes under the global Lipschitz condition were studied by Tian and 
Burrage |22|, Wang et al. [23|, Alcock & Burrage [24] . The tamed Milstein scheme was investigated 
in |25| for SDEs with non-global Lipschitz coefficients. Recently, Higham et al. [26] proposed a new 
Milstein scheme, named as {9, cr)-Milstein scheme, and investigated its strong convergence when 
SDEs hold polynomial growth for the diffusion term. However, the strong convergence rates of 
the implicit Milstein schemes for SDEs with non-global Lipschitz coefficients has not yet been well 
investigated. 
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Stability of numerical solutions is another central problem for numerical analysis. The mean- 
square stability of numerical methods for linear stochastic differential equations have been studied 
by [271 - 129] . For nonlinear SDEs, Higham et al. [30] showed that the BEM and the split-step BEM can 
reproduce the mean-square exponential stability of the exact solution. Recently, without the global 
Lipschitz, Huang m and authors [3211331 presented conditions under which the stochastic theta 
method and split-step theta method can not only reproduce the exponential mean-square stability 
of the exact solution, but also preserve the bound of the Lyapunov exponent for sufficient small 
stepsize, which measures the decay rate of the numerical solutions. However, there is little work 
on the mean-square stability of the theta-Milstein schemes for SDEs with non-Lipschitz continuous 
coefficients. 

The main aim is to examine the boundedness and the convergence rate as well as the expo¬ 
nential mean-square stability of theta Milstein schemes for SDEs with non-Lipschitz continuous 
coefficients. This paper is organized as follows. The next section presents some necessary nota¬ 
tions and preliminaries, and then introduces the stochastic theta-Milstein scheme and the split-step 
theta-Milstein scheme. Section 3 establishes the uniform boundedness of the pth moments of the 
theta-Milstein schemes and shows that for 9 G [1/2,1], these two classes of theta-Milstein schemes 
are bounded in the sense of moment, but for 9 G [0,1/2], the linear growth condition on drift is 
added to obtain the moment boundedness of the theta-Milstein schemes. Section 4 proves that the 
theta-Milstein schemes strongly converge to the exact solution with the order 1. The final section 
investigates the exponential mean-square stability of these two classes of theta-Milstein schemes. 

2 Notations and preliminaries 

Throughout this paper, unless otherwise specihed, we use the following notations. Let | • | denote 
both the Euclidean norm in M”’. If x is a vector, its transpose is denoted by and the inner 
product is denoted by (x,y) = x^y for x,y G M”. a V 6 represents max{a,6} and a A 6 denotes 
min{a,6}. N+ represents the positive integer set, namely, N+ = {1,2,3,...}. Let (D,5,P) be a 
complete probability space with a filtration {5^t}t>o satisfying the usual conditions, that is, it is 
right continuous and increasing while 'So contains all P-null sets. Let be a one-dimensional 

Brownian motion dehned on this probability space. 

Let f,g : M"' i—)• R"' be Borel measurable functions. Let us consider the n-dimensional SDE of 
the form 

dx{t) = f{x{t))dt + g{x{t))dw{t), t > 0 (2.1) 

with initial data x(0) = xq G R”. Assume that / and g satisfy the following assumption: 

Assumption 2.1. Assume that the functions f,gGC^ and there exist constants /r G R and c > 0 
such that for any x, y G R” 

(x - y, /(x) - f{y)) < n\x - y|^, 

Idix) - g{y)\‘^ < c|x-y|2. 
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( 2 . 2 ) 

(2.3) 


Remark 2.1. From conditions (|2.2I1 and (|2.3I) . we have the following monotone condition 


{x, f{x)) V \g{x)\'^ < a + /3|xp, x G M”, 


(2.4) 


where 

« := ^l/(0)P V 2|sr(0)p and /3 := (^ + i) V 2c. 

It is well-known that existence and uniqueness of the global solution of the SDE (mi can be 
guaranteed by the local Lipschitz condition and the monotone condition (or general LV condition) 
(see [M])- Since f,gGC^ implies that / and g satisfy the local Lipschitz condition, the monotone 
condition (|2.4I) can guarantee the existence and uniqueness of the gloal solution for the SDE (12.11) . 
This implies that Assumption 12.II can guarantee the existence and uniqueness of the global solution 
of the SDE (12.11) . 


According to the monotone condition (12.41) . we can present the following lemma (see Lemma 
3.2 in [n]). 


Lemma 2.2. Under Assumvtion \2. 1\ for each p>2 there is C = C{p,T) > 0 such that 

E sup |x(t)|P < (7(1-b E|xo|^). 

Lo<t<T 


In what follows, for the purpose of simplicity, let C represent a generic positive constant inde¬ 
pendent of the stepsize A, whose value may change with each appearance. 

Eixed any time T > 0 and given a stepsize A = T /N for certain integer N, we introduce the 
split-step theta-Milstein (SSTM) approximation 

f yk = Zk + Of{yk)A, ^ 

\ Zk+i = Zk +f{yk)^ + g{yk)^u}k + ^L^g{yk){\A:^Wk\‘^ - A), k = 0,1,2,... ,N. 

where yo = x(0), zq = yo - 6f{yo)A, 9 G [0, 1], = g{x)^ and Awk := w{{k -b 1)A) - w{kA) be 

the Brownian increment. This scheme {zk}k>o can be considered as an extension from the split-step 
theta method developed by our previous work m- It is interesting that the approximation {yk}k>o 
in (12.51) is in fact the stochastic theta-Milstein (STM) approximation 

yk+i = yk + Of{yk+i)A -b (1 - 6)f{yk)A + g{yk)Awk -b g{yk){\Awk\^ - A), (2.6) 

which can be proved by substituting Zk = yk — (^f{yk)A into the second equation in (|2.5p . The 
STM (j2.6p . investigated in [35] for linear scalar SDEs, can be consider as the special case a = 0 of 
the (0, (T)-Milstein scheme developed by Higham et ah [2^. When 9 = 0, {zk}k>o and {yk}k>o are 
the classical Milstein approximation proposed in Milstein [2|. When 9 = 1, {yk}k>o becomes the 
drift-implicit Milstein scheme, which was investigated in |36j . Since theta-Milstein schemes (12.51) 
and (|2.6p are semi-implicit when 9 G (0,1], to guarantee that they are well defined, we restrict the 
stepsize A satisfying 0//A < 1. This, together with the one-sided Lipschitz condition (|2.2p . the 
equation 

y = z + 9Af{y) 


has unique solution y = F^^g{z) for any 2 : G M". 
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3 Uniform boundedness of pth moments 


In order to investigate the boundedness of pth moments and convergence of the two classes theta- 
Milstein approximations, we need the following assumption, which is standard for the classical 
Milstein scheme. 


Assumption 3.1. Assume that the functions f,g £ C'^ and there exists a constant a such that for 
any x, y G M" 

\L^g{x) - L^g{y)\‘^ <a\x-y\^. (3.1) 

Then let us investigate the pth moment boundedness of the theta-Milstein approximations for 
9 G [1/2,1] and 6 G [0,1/2), respectively. 

Theorem 3.1. Let Assumvtions \2. 1\ and \3.1\ hold. 6 G [1/2,1] and let A < A* = l/{26f3). Then 
for each p >2 


and 


E 


E 


sup \Zk\^ 
feAe[o,r] 


sup \yk\^ 
fcAe[o,r] 


< C 


< c. 


(3.2) 

(3.3) 


Proof. By the second equation in (|2.5p . 

kfc+iP = kfeP + l/(yfc)pA2 + i|L^p(pfc)(|Au;fcp - A)|^ 

+2AzJ/(yfc) + {g{yk)^Wk, g{yk){\Awk\‘^ - A)) 

+‘^{zk + f{yk)^,9{yk)^Wk + ^L^9{yk){\^Wk\‘^ - A)) + \g{yk)f\Awk\‘^ 
= \zk\‘^ + 2Aylf{yk) + 2\g{yk)\^\Awk\^ + (1 “ ‘^&)\f{yk)\^Af 
+\\L^9{yk)\^{\Awk\^ - Af + ^(yfc - (1 - 9)zk,giyk)Awk) 

+^{yk - (1 - 0)zk, L^g{yk)i\Awk\‘^ - A)), 

where we used the equation Zk = yk — 9f{yk)A. Note that 9 G [1/2,1]. By (|2.4I) . we have 
kfc+il^ < \zk\^ + 2(a + /3|pfcP)A + 2\g{yk)\^\Awk\^ + ]^\L^9{yk)?‘{\Awk\^ - Af 

2 'Y _ Q Y 

+-Q{yk,9{yk)Awk) - 2—^{zk,9{yk)Awk) + -{yk,L^g{yk){\Awk\^ - A)) 
-^-j^{zk,L^9{yk){\Awk\'^ - A)), 
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which implies 


k k 

kfc+iP < \zo\‘^+ 2aT + f3Aj2\yi\‘^+ ‘^'^\9{yi)\‘^\^Wi\^ 

i=0 i=0 

1 fc 2 ^ 

+Q'^{yi,L^9iyi)i\^Wi\‘^ - ^)) + g'^{yi,9{yi)^wi) 

i=0 i=0 

k k 

+ '^'^\L^g{yi)\^{\Awi\‘^ - Af 

i=0 i=0 

k 

-—0-^^^i^^^9{yi){\^Wi\^ - ^))- 
i=0 

Recall the elementary inequality: for xi, ■ ■ ■ , x/ > 0, p > 1, / = 1, 2,..., A^, 




2 = 1 


2=1 


We therefore have 


8P 


^ k k 

< (|^oP + 2aTr + /3W(^|y,|2)%2P(^|5(yi)P|Au;i| 

2=0 2=0 
k k 

+2^1 '^{yi,L^9iyi)i\^'Wi\‘^ - A)) ^ + 4^^! 'Y^{yi,g{yi)Awi 

i=0 
k 

+ 2^1 '^{zi,g{yi)Aw 


i=0 


+ '^{zi,L^g{yi){\Awi\^ - A)) 


i=0 

k 

+ |L^p(yi)p(|Au;i|2 - Af 

i=0 

p 


(3.4) 


1=0 


Note that yi is Jj- measurable while Awi is independent of Hence, for any integer m < N, 


k 

Ie( sup V 

^0<k<m~^ 


\g{yi)\^\Awi\^Y 




sup '^\g{yi)Awi\^P 

<k<m • n 
- — 2 = 0 

n 

= NP-^E[Y,\9{yi)^m\^^' 

2 = 0 
m 

= NP-^J2^\g{yi)\^PE\Awi\^P 
2=0 
m 

< CAj]E[« + /3|yil2]P 

2=0 

m 

< C'A^E[aP + /3^'|2/,|2?>] 

2=0 

m 

< c + CAj2n\yi\‘'^], 

i=0 


(3.5) 
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where we also used (12.4p . Using the Burkholder-Davis-Gundy inequality, we have 


E 


K m 

sup I - A)) n < CE\^\yi\^\L^g{yi)\ 


*- 0</c<m 


,2^2 


lP/2 


and 


i=0 


i=0 


< CAP(m + J;E[|2 /,P|l15(0)P + a\y,\^f/^ 

z=0 

m 

< C'A^’/^ + C'AP/^+^ E (3.6) 

i=0 


E 


sup E yTg{yi)^Wi 


*- 0<k<m ' ^_Q 


< 


m 

CE[Y^\y,\^\g{y,f 


A 


i=0 


p/2 


< CAp/2(^ +i)p/2-1e E |yir[a + /3|y,|V 

m 

< C + UA^[l+E[|y,l2p], 

i=0 


E 


sup 

0<fc<m ' ^_Q 


g{yi)Awi 


< 


p/2 


- C'IE[^|2:ip|fi((yi)pA 

m 

< CAP/2(m + l)P/2-iE Ei^<n“ + /3|yilV 

i=0 

m m 

< CAY,[l+E[\y,\^P] + CAY,n\z^\^n, 


i=0 


i=0 


E 


K m 

sup \^zf L^g{yi){\Awi\^ - A))^] < CE\^\zi\^\L^ g{yi)\ 


0<A:<m ^_Q 


,2a2 


p/2 


i=0 


(3.7) 


(3.8) 


< UA?’(m + l)P/ 2 -i^E[|^^| 2 |^i^(o)| 2 ^^|y^| 2 ]p /2 

i=Q 

m m 

< U + CA^EOy.pPJ + CA J]E[|zi|2p]. (3.9) 


i=0 


i=0 
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Similarly, we have 


sup V |L^ 5 (yi)| 2 (|Au;ip - A)2] < sup ^ |L^ 5 r(yi)(|Au;*^ - A)|^p 

0<fc<m ^_Q ^ ^0<fc<m^_Q 

m 

= NP-^E[Y^\L^g{y,){\Aw,\^-A)\^P 

i=0 

m 

i=0 

m 

= cAY,nL^gm^+<y\yr?r 
2 = 0 
m 

< CAY,E[\L^gm^P + a\y,\^P\ 

2 = 0 

m 

< c + CAY,nyi?^]- 

Combining (|3.5p - (l3.10l) with (13.41) yields 


i=0 


E 


sup \zk\^P < C + CAVE[|y,|2P] + CAVE[|y,l2P] 
0<fc<m+l J 


i=0 

m 


i=0 


< c + 


i=0 


i=0 


0<k<i 


Using Zk = yk- Gf{yk)^ and (12.41) . we have 


kfcl^ = lykl"^-2eylf{yk)A + 9‘^A^\f{yk)\‘^ 

> (l-20A/3)|yfc|2-20Aa, 


(3.10) 


CAVEf sup |yfc|H +C'AVe[ sup (3.11) 

^ ^0<k<i J \-n<h<i J 


that is, 


(l-20A/3)|yfc|2 < \zk\^ + 2eAa, 
which together with (13.lip implies that for A < A* 


E 


sup \Zk 
'■ 0<fc<m+l 


|2p 


<c + 


m 

CA^Ef sup {zkl"^^ 
^ ^0<k<i 


i=0 


Using the discrete-type Gronwall inequality and noting that [m -|- 1)A < T give 


E 


sup \Zk\^^ 
0<fc<m+l 


< c. 


This together with (|3.12l) gives the desired assertion (13.31) . 


(3.12) 


(3.13) 

□ 


The following theorem gives the moment boundedness for 0 € [0,1/2). 




















Theorem 3.2. Let Assumptions \2. 1\ and \3.1\ hold. 0 G [0,1/2) and let A < Ai = l/(20/3) (A 
if 9 = 0). If function f satisfies the linear growth condition 


then for each p>2 


and 


Proof. By (|2.5I) . 


\f{x)f < K{1 + |xp), 


E 


sup \yk\^ 
fcAe[0,T] 


< C 


E 


sup 

fcAe[o,T] 



< C. 


\zk+i\^ = \zkf+ 2{zk,f{yk)t2^ +9{yk)^Wk +-L'^g{yk){\Awk\^ - A)) 


+ 


/(yfc)A + g{yk)Awk + ^L^g{yk){\Awk\'^ - A) 


Note that 


Zk = yk- 9f{yk)A. 


Substituting this into (I3.17P produces 


kfe+ip < \zk\^ + 2{ykJ{yk)Ap g{yk)Awk + ]^L^g{yk){\Awk\^ - /A)) 

+0'^^^\f{yk)? + 2|/(yfc)A + g{yk)Awk + g{yk){\Awk\^ - A)p 

< \zk\^ + 2ylf{yk)A + 2ylg{yk)Awk + ylg{yk){\Awk\^ - A)) 
+7AV(2/fc)P + 6|fif(2/fc)Awfcp + 3|L^fi((?/fc)(|Auifc|^ - A)p 

< jzfcp + (2a + 7|/(0)|2A)A + (2/3 + IKA)\yk\'^A 
+‘^yk9{Vk)^Wk + ylL^g{yk){\Awk? - A)) 

+Q\g{yk)Awk\^ + ?,\L^9{yk){\Awk\^ - A)p, 


which implies 


k 

kfc+iP < koP + (2a + 7|/(0)|2A)r + (2/3 + 7ii:A)A^|y,|2 

2=0 

k k 

2=0 2=0 

k k 

+6^ b(2/i)Au;ip + g{yi){\Awi\'^ - A)p. 

2=0 2=0 


L = OC 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 


(3.19) 
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Similar to (13.41) . we have 


6P 


^-\zk+i\^^ < (|zoP + (2a + 7|/(0)|2A)Tf+ (2/3 + 7KA)W(^|y. 

i=0 

k k 

+2^1 ^ + I 

i =0 2=0 

k k 


i=0 

Combining (l3.5D - (l3.7p . (I3.1UI) with (|3.2UI) yields 

E 

which together with (I3.12p implies 


i=0 


sup \zk\'^P < C + CA VE[|yipP], 
'-0<k<m+l J 


i=0 


E 


sup \Zk 
0<fc<m+l 


|2p 


<C + CA^E[|yil2P] < C + CA^E[ sup \zk\‘^P] 


i=0 


i=0 


0<k<i 


(3.20) 


Note that (m + 1)A < T. Using the discrete-type Gronwall inequality gives the desired assertion 
(I3.16p . Similarly, using (|3.12l) and (|3.16l) gives another desired assertion (|3.15l) . □ 


Remark 3.1. For 6 £ [0,1/2), the linear growth condition p3.14p may be necessary for guaranteeing 
the moment boundedness of the theta-Milstein schemes. This fact can be found in [321 Lemma 3.2] 
and [Silo]. 


4 Convergence of the theta-Milstein approximations 


Let us now introduce appropriate continuous-time interpolations corresponding to the discrete 
numerical approximations. More accurately, let us define the continuous solutions z{t) and y{t) for 
t G [tk,tk+i) as follows 

z{t) = z{tk) + {t- tk)f{yk) + g{yk){w{t) - w{tk)) + ^L^g{yk){\w{t) - w{tk)\^ - {t - tk)), 
y{t) = FA,e{z{t)), 

(4.1) 

with 2;(0) = zq = yo — 0A/(yo), where tk = k/S.. Note that z{t) = y{t) — 9f{y{t))A. Then the 
continuous-time approximations z{t) and y{t) are measurable. In our analysis it will be more 
natural to work with the equivalent definition of z{t) 


f 


z{t) = z(0) + / f{y{s))ds + / g{y{s))dw{s) -h 



L^giy{u))dw{u)dw{s), (4.2) 


where s = for s G [tk,tk+i)- Note that z{tk) = Zk and y{tk) = yk- We refer to z{t) and y{t) as 
the continuous-time extensions of the discrete approximations Zk and yk, respectively. 
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In what follows, we also use Taylor’s formula frequently. If a function h : M” —)• M” is twice 
differentiable, the following Taylor’s formula holds 

h{z{s)) - h{z{s)) = h'{z{s)){z{s) - 2 ( 5 )) + Rs{h), (4.3) 

where Rs{h) is the remainder term dehned by 

Rs{h) = [ {1 - r)h”{z{s)+ r{z{s) - z{s))){z{s) - z{s),z{s) - z{s))dr. (4.4) 

Jo 

Here for any tt, u G M” the derivatives have the following expression 

'•'(■xx) = E '‘"(-xx,*) = E 

2=1 ^J = l 

Replacing z(s) — z{s) in ()4.3p with ()4.ip and rearranging lead to 

h{z{s)) - h{z{s)) = h'{z{s)) g{y{s))dw{s)J + Rs{h), (4.5) 

where 

Rs{h) = h'{z{s)) (^(s - s)f{y{s)) + ^L^g{y(s)){\wis) - 'u;(s)p - (s - s))^ + Rsih). 

In order to obtain the strong convergent rate, we also need the following assumption. 
Assumption 4.1. There exist positive constants D and q such that for all x G MT 


\f{x)\y\r{x)\<D{i + \xn 


(4.6) 


and 


\9"ix)\ < D. 


(4.7) 


Let us hrstly give the convergence theorem of the theta-Milstein schemes for 6 G [1/2,1] as 
below. 


Theorem 4.1. Let Assumptions l2Al 13.11 and \4-l\ hold, 6 G (1/2,1] and let A < A* 
Then for any p>2, 


E 


sup \x{t) 
te[0,T] 


-zit)f 


< CAP 


and 


E 


sup \x{t)-y{t)\P 
te[o,r] 


< CAP. 


To prove Theorem 14.11 we need the following lemmas. 


l/(20/3). 

(4.8) 

(4.9) 


Lemma 4.2. Let the conditions in Theorem\4.1\ hold. Then for any p > 2, 


sup 

0<k<N 


nfiyk)\^ 


V 


sup 

0<fc<Ar 


nnyk)\ 


p 


V 


sup 

0<k<N 


nfiykW 


< c 


(4.10) 


and 


sup 

0<k<N 




V 


sup 

0<A:<Ar 


nL^9{yk)\ 


p 


< c. 


Moreover, replacing yk by Zk, (I4.10p and (|4.11l) are also true. 


(4.11) 
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Proof. By the polynomial growth condition (14.611 . the global Lipschitz condition for g and L^g and 
Theorem 13.11 give the desired assertions immediately. □ 

Similarly, we can use the polynomial growth condition (14.6|) . the global Lipschitz condition for 
g and L^g and Lemma 12.21 to obtain the following lemma, whose proof is omitted. 


Lemma 4.3. Let conditions in Theorem\4^. 1\ hold. Then for any p >2, 


sup 
0<t<T L 


E\f{x{tW 


V sup 
0<t<T L 


nnx{t)w 


V sup 
0<t<T 


E\f{x{tW 


< C 


(4.12) 


and 

sup []E|ff(y(t))|H V sup \E\L^g{y{t))\P 

0<t<T L J 0<t<T L 

Lemma 4.4. Let conditions in Theorem\4.1\ hold. For any p >2, 


< C. 


E 


sup \z{t)\P 
'-te[o,r] 


VE 


sup \y{t)\P 
4e[o,T] 


< C 


and 


E 


sup \z{t)-y{t)\P 
'-te[o,T] 


< CAP. 


(4.13) 


(4.14) 


(4.15) 


Proof. For any p>2 and t G [0, T], by (14.2p . 


E(" sup |z(s)|^') < 4 P ^E|z(0)|^ + 4^ ^E f \f{y{s))\ds +4^ ^E^ sup f g{y{s))dw{s) 

s£[0,t] ' -In ^ ocTn fl .In 


se[o,t] Jo 


+4P-^e( sup 
^se[o,t] 


L^ g{y{u))dw{u)dw{v] 


0 Jv 


(4.16) 


Using Lemma 14.21 we have 


E 


/ |/(y(s))|ds < CE \fiy{sWds<C. 

/o Jo 


(4.17) 


Applying the Burkholder-Davis-Gundy inequality and Lemma 14.21 yield 

ri - p/2 


e( sup [ g{y{s))dw{s) H < CpE( f |c/(y(s))pds)^ 
^se[o,t] Jo ' ^Jo ' 

< C \\\g[y{s))fds<C 
Jo 


(4.18) 


and 


E( sup 

' sG[0,t] 



0 J V 


L^ g{y{u))dw{u)dw{v)Y < CpE(^ J J L^g{y{u))dw{u) 

CAp!^ [ E\L^g{y{v))\Pdv 
Jo 


2 ^p/2 

dv 


< 

< C, 


(4.19) 
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where Cp is a constant dependent on p. Therefore, (I4.17I) - (I4.19I) and (I4.16P produce 


which together with (j3.12l) gives 


E 


E 


sup | 2 :(s)l^ < C, 
se[o,T] -I 


sup |i/(s)|P 
se[o,T] 


< C. 


Note that z{t) = y{t) — 9f{y{t))A. Hence (I4.15P follows from (I4.14p . 

Lemma 4.5. Let conditions in Theorem \4- 1\ hold. For any p>2 and t S [0,r], 

E|x(t) -x(t)r < CAP/^ 

and 

E\z{t) - z{i)\P < CAP/\ 

Proof. By ()2.1I) and the Burkholder-Davis-Gundy inequality, we have 

E|x(t) — x(t)|^ = CEI f f{x{s))ds\P + CE\ f g{x{s))dw{s)\P 

Ji Jt 


< CAP- 


p/2 

'W E\f{xis))\Pds + CE / \g{x{s))\‘^ds < 

Jt Jt 


where we also used Lemma I4..S1 By (|4.ip . we have 

E\zit) - zii)\P = CAPE\f{y{iW + CE\g{yii)){w{t)-wiiW 

+CE\L^g{y{i)){\w{t) - w{i)f - {t - i))\P 
< CAPE\f{y{i))\P + CAP/^E\g{y{i)W + CAPE\L^giym^ ■ 

Then we obtain (j4.2ip from Lemma 14.21 

Lemma 4.6. For any p>2 and t G [0,T], 

E\Rt{f)\P < CAP. 

Proof. By the definition of Rt{f), we have 


□ 

(4.20) 

(4.21) 


□ 


(4.22) 


E\Rt{f)\p < AP3P-^E\f{zmmt)w+^^-^mt{fw 

+AP3P-^^E\f{zmL^giyii)W^\wit) - wii)\^ - (t - i)) 

By Lemmas 14.41 and 14.51 as well as Assumption 14.11 we can obtain E\Rt{f))\P < CAP. By Lemma 
14.21 we can obtain the desired assertion (14.221) . □ 
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Proof of Theorem \4-l\ Let e{t) = x{t) — z{t). By (|4.2I) and (12.11) . we have 


pi pi p s 

e(t) = e(0)+ / [f{x{s))-f{y{s))]ds+ / [g{x{s)) - g{y{s)) - / g{y{u))dw{u)]dw{s). (4.23) 

Jo Jo Js 


For the purpose of simplification, define G{s) = g{x{s)) — g{y{s)) — Jf L^g{y{u))dw{u). Applying 
Ito’s formula to |e(t)p gives 

|e(t)|2 = |e(0)p+/' \G{s)\‘^ds + 2 f {e{s),G{s))dw{s) + 2 f {e{s), f{x{s)) - f{y{s)))ds 

Jo Jo Jo 

< 9‘^A‘^\f{x{0))\‘^ + f \G{s)\^ds + 2 f {e{s),G{s))dw{s) + 2g f \e{s)\‘^ds 

Jo Jo Jo 

+2 [ {e{s),f{z{s))-f{y{s)))ds + 2 [ {e{s), f{z{s)) - f{z{s)))ds 
Jo Jo 

= : 0^A^|/(x(O))|^ + J\it) + ^ 2 (f) + 2/i f |e(s)|^(is + Jsit) + (4.24) 

40 

Applying Ito’s formula for g{x{s)) gives 

g{x{s)) = g{x{s)) + ( g{x{u))dw{u) + [ (g o f){x{u)))du, (4.25) 

J S J S 

where {g o f){x{u))) = g'{x{u))f{x{u)) + ^trace{g{x{u))'^g''{x{u))g{x{u))}. Substituting (j4.25p 
into G{s) and using the Burkholder-Davis-Gundy inequality and Holder’s inequality yield that for 
any p > 2, 

IE|G(s)l^’ < AP-^E\gix{s))-g{y{sW + 4P-^E [\L^gix{u)) - L^gix{u))]dw{u) 


+ 4 P-IE 


< GE|x(.s) — y(s)|^ + CE j \x{u) — x{u)^]ds 

p/2 


[L^g{x{u)) - L^g{y{u))\dw{u)^ + 4P ^E j {g o f){x{u)))d 

. ,0. . P/2 


+CE 


J \x{u) — y{u)\‘^ds 


+ CAP 


< GE|a;(.s) — 2 ;(s)|^ + GE| 2 :(s) — y(s)|^ + ^E / \x{u) — x{u)\^ds 

J S 

PS PS 

+GAP/2-1E / \x{u)-z{u)\Pds + GAP/^-^E \z{u) - y{u)\Pds + GAP 
J S J S 

< GE|e(s)|P + GAP, 

where we also used Assumptions 12.11 and 13.11 as well as Lemmas 12.2114.2114.51 Hence, 


E 


sup |Ji(t)|P/2j = GE / |G(s)|Pds < G / E|e(s)|Pds + GA^. (4.26) 

te[0,r] -I Jo Jo 


Using the Burkholder-Davis-Gundy inequality and the Holder inequality yields that for any p > 2 
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and r G [0, T] 


E 


sup |J 2 (t)|^^^ = 2 ^/^E sup f {e{s),G{s))dw{s) 

tS[0,r] ts[0,r] Jo 


p/2 


p/4 


tS[0,r] Jo 

< 2P/‘^CpE(^ j \{e{s),G{s))\^ds 

pT 

< 2 ^/^CpEf sup le(s)|^ / \G{s)\‘^ds 

^ sG[0,r] Jo 


p/4 


where Cp is a constant dependent on p. Recalling the fundamental inequality; 2ab < kio^ + ^ 6 ^ 
for any a, b,Ki > 0 yields 

E 




2Pcl 

< AtiE 

sup e(s)|^ 

+ 


^sS[ 0 ,r] 

Kl 

VI 

sup |e(s)|^ 

+ c [ 


^se[ 0 ,r-] 

Jo 


sup \J2{tW^‘^ 

4e[0,r] 


Similarly, using Holder’s inequality, Assumption 14.11 and Lemma 14.41 we have 


vP/2 


(4.27) 


E 


sup |J 3 (t)|^/^ < GE f \e{s)\Pds + GE f \f{z{s)) - f{y{s))\Pds 

'-tG[0,r] ^ Jo Jo 

pr pr 

<G E|e(s)|^ds + C / ■\/E(1 + | 2 ;(s)|'? + |?/(.s)|'J)2pE|z(s) — y{s)\'^Pds 

Jo Jo 

< G [ E\e{s)\Pds + GAP. (4.28) 

Jo 

Now, let us estimate J 4 . Using (|4.5I) produces 

Mt) = {e{s),f{z{s))(^J^g{y{s))dw{ufj+Rs{f)'jds 

‘J {e{s),f{z{s))(^j^g{y{s))dw{ufj'jds + 2j {x{s) - z{s), Rs{f))ds 

J 4 i{t)+ f\e{s)fds+ f\Rs{f)\^ds, (4.29) 

Jo Jo 

where J 4 i{t) = 2 jQ{e{s), f^{z{s)){JJ g{y{s))dw{u)))ds. It is easy to deduce from Lemma HR] that 


= 2 


< 


E 


sup 

'-te[0,r] 


f\Rs{f)\^ds 

Jo 


p/2l 


< GAP. 


(4.30) 


Noting that e(s) = e(s) + J~[f{x{u)) — f{y{s)]du + G{u)dw{u), we have 

J 4 i(t) = 2 f f {e{s), f {z{s)){g{y{s))))dw{u)ds 

Jo Js 


0 Js 

+2 j G{u)dw{u), f {z{s))(^ j ^ g{y{s))dw{u)^'^ds 

R2 ' ^ ^ 


'lo^L 9{y{s))dw{v)j'^ds 

= : J4n{t) + J412{t) + J413{t)- 
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For any t G [0,T], let kt = max{A; > 0,tk < t}. Define s := ti for ti-i < s < U and s := t for 
tkt < s < t. Applying Fubini’s theorem, we can obtain 


kt pti ps 

J4ii{t)/2 = ^ / / {e{ti-i),f'iz{ti_i)){g{y{ti-i))))dw{u)ds 

+ [ [ {e{tkt)J'{z{tk,)){g{y{tkt))))dw{u)ds 

= [ f {^i^i-i)J'i^iti-i)){g{y{ti-i))))dsdw{u) 

+ [ [ {e{tkt), f {z{tkt)){ 9 {y{tkt))))dsdw{u) 

Jtkt du 

kt nti 

= '^ (ti-u){e{ti_i),f'{z{ti-i)){g{y{ti_i))))dsdw{u) 

i=l JU-i 

+ [ (t-u){e{tk,)J'(z{tk^))(g{y(tkt))))dsdw{u) 

dtkt 

= [ (u-u){e{u),f{z{u)){g{y{u))))dw{u). 

Jo 

For any K 2 > 0, using the Burkholder-Davis-Gundy inequality, Holder’s inequality and Lemma 14.21 
gives 

,P/4 


E 


'-te[0,r] 


sup \JAiiit)\P^‘^ < CpE( / 2\{u - u){e{u), f{z{u)){g{y{u))))\‘^du 


0 


< ‘2\e{u)\‘^A‘^\f{z{u)){g{y{u)))\‘^du 


p/4 


< CpEh sup le(s)|^A^ / \f'{z{u)){g{y{u)))\'^du 
^ se[o,r] Jo 


p/4 




cl2P/^ / 

< K 2 E 

sup e(s)P 

+ ^ APE( 


^sS[0,r] 

K2 V 

VI 

sup e(s)|^ 

+ CAP, 


^sS[0,r] 



p/2 


(4.31) 
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Similarly, 


E 


sup < CE [ / [ G{u)dw{u),f{z{s))( [ g{y{s))dw{u))) 

te[o,r] J Jo ^Js ^Js 


p/2 


ds 


< C WE 


< C WE 


E 

ps 

/ G{u)dw{u 
J S 

p 

E 

y 

f'izis)){ 9 {y{s)))dw{u) 

s 

^ds 

E 

ps 

/ G(m) 2(iu 

J S 

p/2 

E 

J 

[ \r{z{s)){g{y{s)))\^du 

S 

p/2 

ds 


< c ^/AP-2 / E\Giu)\Pdu I E\r{zis)){giyim\Pduds 


< 


< 


rr 

C / VAP[E|e(s)|P + AP]ds 
Jo 

G [ E\e{s)\Pds + GAP. 

Jo 


(4.32) 


In order to estimate J 413 , we divide it into the following three parts 

r*t , rs 


Jiloit) = 2j^(^J^[f{x{u))-f{x{u))]duJ'{z{s))(^J_g{y{s))dw{vfj'^ds 

^‘^JoiJ. 9{y{s))dw{v)'j'^ds 

lo^L 9 {y{s))dw{v)j'^ds 


+2 

/o ' Js 
= : Ii{t) + l2{t) + Is{t). 


(4.33) 


Using the Burkholder-Davis-Gundy inequality and Holder’s inequality give that for any p > 2, 
e[ sup 

■te[0,r] 

p/2 


< 


CE J [f{x{u)) - f{x{u))]du, f {z{s))(^ g{y{s))dw{u)j'j 


ds 



<G WE 

'0 


[f{x{u)) - f{x{u))]du E / f{z{s)){g{y{s)))dw{u) 


ds 


< G / WAP-i 


mi^iu)) - f{x{u))\PduAP/‘^-^ E\f'{z{s)){g{y{s)))\Pduds 

pr _ 

<G Vap- 1 aAp/2Ap/2- 1 Ads < GAP, (4.34) 

Jo 
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where we used Lemma 14.21 and 14.51 Similarly, 


E 


'-te[0,r] 



p/2 


ds 


< 


< 


C ^APE\f{z{s)) - /(y(s))|PAP/2-i ^ E\f{z{rs)){g{y{s)))\Pduds 

pr _ 

C / VAPAPAP^ds 
Jo 


< CAP. 


(4.35) 


Applying Fubini’s theorem obtains 

rt rs 


h{t) = 2 / / {s - s){f{x{s)) - f{z{s))J'{z{s)){g{y{s))))dw{v)ds 

Jo Js 

[ [ (s- s){f{x{s)) - f{z{s)), f'{z{s)){g{y{s))))dsdw{v). 
Jo Jv 


= 2 


Therefore 

E 


sup \h{t)\P/^ 

'■te[0,r] 


< CE 


< 


< CAP 

< CAP. 

By (I4.34I) - (I4.36|) . we obtain from (I4.33D 

E 


I (s - v){f{x{v)) - f{z{v)), f {z{v)){g{y{v))))ds 

V 

[ W{x{v)) - f{z{v)), f'{z{v)){g{y{v)))) 

Jo 


2 ^p/4 

dv , 


2 ^p/4 

dv 


p/2 


(4.36) 


sup \JAiz{t)\P/‘^ < CAP, 

te[0,r] J 


which together with (14.311) and (I4.32p implies 

E sup |J4i(f)|^/^ < 2)P/‘^~^(^ sup |J4ii(f)|^/^ 

tS[0,r] ^ ^tS[0,r] 


+ E 


sup \JAl 2 {t)\P^‘^ 
te[0,r] 


+E 

Combining (I4.37|) . (I4.30p with (I4.29P gives 


sup |J4i3(i)r^^ 
'-te[0,r] 


sup |e(s)|^ 

'-se[0,r] 


+ C 


l\\e{s)\Pds 

Jo 


+ CAP. (4.37) 


E 


sup |J 4 (f)|^/^ 

te[0,r] 


< 3P-^K2'E 


sup |e(s)|P +C E\eis)\Pds + CAP + C / E\e{s)\Pds. 
se[0,r] -I Jo Jo 
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Let K = (ki + 3^ ^. Combining (I4.26p . (14.271) . (I4.28p . (I4.38P and (|4.24l) produces 


E 


sup le(t)|^ 
te[o,r] 


< 


6P/2-1Ve[ sup \Jl{t)\P/^ 
'-te[o,r] 


+ CAP + C 


1=1 


f le(s)|P(is 
Jo 


< kE 

sup 6 ( 5 )/^ 

+ CAP + C 

pr 

/ E 

sup \e{u)\P 


■sS[0,r] 


/o 



ds. 


Therefore, letting ki and K 2 be sufficiently small such that k < 1, we have 

ds + CAP. 


E 

sup \e{t)\P 

< c 

pr 

/ E 

sup |e(s)|^ 


^t£[0,r] 


/o 

■uS[0,s] 


(4.38) 


Thus, the result ()4.8I) follows from the Gronwall inequality. The desired assertion (j4.9p follows from 
and Lemma 14.41 □ 


For 9 E [0,1/2], we can use the similar techniques to obtain the following convergence theorem. 
Its proof is omitted. 


Theorem 4.7. Let Assumptions \3.1\ and \4-l\ hold, 6 G [0,1/2) and let A < Ai = l/(20/3) 

(Ai = oo, if 9 = 0). If the function f satisfies the linear growth condition (I3.14p . then for any 
p >2 we have 


E 


sup \x{t) 

t€[0,T] 


zit)f 


< CAP 


and 


E 


sup \x{t)-y{t)\P 
te[o,T] 


< CAP. 


Remark 4.1. Convergence for 9 G (1/2,1] in Theorem 14.11 are partially based on the Assumption 
O This does not imply that condition (12.31) in Assumption 12 .1 1 is necessary. For example, Kloeden 
and Neuenkirch m showed that the semi-implicit Milstein scheme (STM scheme with 9 = 1) 
converges to the exact solution for the Cox-Ingersoll-Ross process, whose diffusion term is non- 
Lipschitz continuous. 


Remark 4.2. Higham and his coauthors mm showed that the split step backward Euler and 
backward Euler schemes strongly converge to the exact solution with the standard order 0.5. In 
fact, by the similar skills used in the proof of Theorem 14.11 we can show that the convergence rate 
0.5 holds for theta-Euler schemes (split-step theta and stochastic theta schemes) with 9 G (1/2,1]. 
Moreover, for the SODE with additive noise, that is, g{x) is a constant, then L^g{x) = 0, the 
theta-Milstein schemes are degenerated to the theta-Euler schemes. Hence, our convergence results 
show that the theta-Euler schemes {9 G (1/2,1]) is convergent with order 1, which is a new result 
for the SODE with one-sided Lipschitz and polynomial Lipschitz continuous drift. 


In this section, although the strong convergence results for the theta-Milstein schemes are 
established under Assumptions 12.11 13.11 and 14.11 the following exponential mean-square stability 
only depends on the one-sided Lipschitz condition (12.2p and the local Lipschitz condition. 
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5 Stability analysis of the theta-Milstein schemes for SDEs 


For the purpose of stability, assume that /(O) = g(0) = 0. This shows that (|2.ip admits a trivial 
solution. Then inequality (13.111 in Assumption 13.II becomes 

\L}g{x)\'^ < (7|xp. (5.1) 

In this section, we assume that / and g satisfy the following local Lipschitz condition, which is 
classical for the nonlinear SDEs. 


Assumption 5.1. (Local Lipschitz condition) f and g satisfy the local Lipschitz condition, that 
is, for each j > 0 there exists a positive constant Kj such that for any x,y £ M” with | 3 :| V \y\ < j, 

\f{x) - f{y)\ V \g{x) - g{y)\ < Kj\x - y\. (5.2) 

To investigate stability of numerical approximations, let us firstly give the stability criterion of 
SDE (12.11) (also see [M]i: 

Theorem 5.1. Let Assumvtion I5.il hold. If there exists a positive constants 7 such that for all 
x,y £ M”, 

2x^f{x) + \g{x)\'^ < - 7 |a:|^, (5.3) 

then the solution of m has the property 

E|x(t)p < C{xo)e~'^^, (5.4) 


where C{xo) is a positive constant dependent on the initial data xq. 
Under the one-sided Lipschitz condition (12.2p . we dehne 


A* 


1 

00 , 


9 £ ( 0 , 1 ), and // > 0 , 
0 = 0 or // < 0 . 


When A < A*, the two classes of theta-Milstein schemes are well defined. Moreover, we have the 
following stability result. 


Theorem 5.2. Let the conditions in Theorem ED hold. Under condition (EH) and the one-sided 
Lipschitz condition ( 1 ^. the following two stability results hold. 

(i) Let 0 G [0,1/2]. If there exists a positive constant K such that for any x £ M"’, 

|/(x)p < iL|x|^, (5.5) 


then for any A < A! = --- A A*, k £ N+, 

■' ^ ^ {1 -2e)K + 0.5a + 


where 7 a 


< C(xo)e-^^'=^ and Ely^j^ < C7(xo)e-^^'=^, 

1, /, 7 -(l-26i)ArA-0.5o-A ^ 

— log 1 - 4-^- L ——-A 

A V (1 + 0A/A)2 
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{ii) Let 6 G (1/2,1]. Then for any A < A 2 = — A A*, k ^ N+, 

a 

E\zkf <C{xo)e-^^’^^ and E\yk\^ < C{xo)e-'^^’^^ 
(26»- l)(7-0.5aA)A \ 


where 7A = log (1 - _ 1 + _ o. 5 aA)A ^2 j 

Proof. For any large I > 0, define the stoping time 

A; = inf{i > 0 : \yi\ > I or \zi\ > 1}. 

It is observed that \yk-i\ < I, \zk-i\ < I ioi k £ [0, A/]. By (j2.5jl . 

kfcP < 4.{\zk-i\‘^ + |/(yfc_i)pA2 + \g{yk-i)\‘^\Awk-i\‘^ + ^\L^g{yk-i)\‘^\\Awk-i\‘^ - Ap^. 


(5.6) 


Then we can obtain that 


E[\zk\%oM{k)] < 4/2 + 4E[|/(yfc_i)|2l[o,A,](A:)]A2 + 4E[|5(yA:-i)|"|Aw;fc_ipl[o,A,](A:)] 
+E[\L^g{yk-i)\‘^\\Awk-i\'^ - Apljo^;^^](A:)] 

1 /2 

< 4/2 + 4E[|/(2/fc_i)|2l[o^;,^](A:)]A2+4(E[|<7(yfc-i)|"l[o,A4(^)]E[|Aa;fc_i|"]) 

+ (E[|Li5(yfc-i)l'l[o,A4(^)]E[|(Aa;fc_i)2-Ar])'^', 

where we used Holder’s inequality. Making use of the property E[|Au;fcp*] = (2i — 1)!!A*, we 
have E[|Ar(;fc_i|^] = 3A^ and E[|(Ar(;fe_i)^ — A|^] < 2^(E[|Arefc_i|®] + A^) = 2^(7!! + 1)A^, where 
(2i — 1)!! = (2z — l)(2z — 3) • • • 3 • 1 for i = 1,2,.... By Assumption 15.11 and (15.ip . it is easy to 
deduce that |/(yfc-i)pl[o,A 4 (A:) < Kff, \g{yk-i)\^l[o^Xi]{k) < KfL^ and \L^g{yk-i)\'^l[o^Xi]{k) < 
a‘^1^. Therefore, 

E[\zk\‘^l[o,Xi]{k)] < 41^+ 4:K^l^A^+ 8Kfl^A+ 4VWTTafA‘^ =: K{1), (5.7) 

which together with Zk = yk — &Af{yk,yk-m) and 2x'^f{x) < — 7 |xp yields 

= |yfc|^ - 29Aylf{yk) + \f{yk)\‘^A‘^e'^ > (1 + 7 ^A)|yfcp. (5.8) 


Hence, 

E[\yk\%oMik)]<K{l). (5.9) 

By 

\zk+i\‘^ = \zk\‘^ + \ f{yk)\‘^A‘^ + ^\L^g{yk){\Awk\‘^ - A)p + {g{yk)Awk,L^g{yk){\Awk\^ - A)) 
+2{zk + f{yk)A,g{yk)Awk + ^L^g{yk){\Awk\^ - A)) + \g{yk)\‘^\Awk\^ + 2Azlf{yk) 
= \zk\^ + 2Aylf{yk) + \g{yk)\^A + (1 - 2e)\f{yk)\^A^ + ]^A^\L^g{yk)\^ + ruk, 
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where mk = 2{yk + (1 - d)f{yk)A,g{yk)Awk + ^L^g{yk){\^Wk\‘^ - A)) + \g{yk)\‘^{\Awk\‘^ - A) + 
{g{yk)Awk, L^g{yk){\Awk\‘^ - A)) + l\L^g{yk)\'^[{\Awk\‘^ - A)^ - 2A\ Using condition ([5A]) and 
the coupled condition (15.3p . we get 

kfc+iP < \zk\^ - lA\yk\^ + (1 - 29)A^\fiyk)\^ + 0.5A^a\yk\^ + (5.10) 

Case (i) 9 G [0,1/2): Note that the linear growth condition of / and (j5.3p imply 

\g{x)\^ < —2x^f{x) < |xp + l/(x)p < (1 + K)\x\‘^. (5-11) 

From conditions (15.ip . (15.51) and (15.lip , there is a positive constant C such that 

\'mk\ < [C + C\Awk\^ + C\Awk\^\\yk?‘■ (5.12) 

Noting that yk and l[o,Ai](^) 'SkA- measurable while Awk is independent of JfcA) then we obtain 
from (|5.9I) that E[mfcljo,Ai] (^)] = 0- The linear growth condition (15.51) gives 

\zk+i\^ < \zk\^ + [(1 - 29)KA + 0.5ctA - 'y]A\yk\^ + mk- (5.13) 

Then for A < A^ and sufficiently large /, by (I5.10j) and using linear growth condition (15.5p . we have 

(fcAAi)-l (fcAAi)-l 

l^fcAAil^ < 12 :o1^ + [(1-26»)A'A + 0.5(tA-7]A ^ \yi\^ + ^ rm 

2=0 2=0 

k-1 

< |2o|^ + X]^[o,Ai](*)"ii- (5-14) 

i =0 

Hence, for A < A^, E[|zfcAAiP] < lE|zoP = E|yo — ^A/(yo)p =: K. Note that 

n\zx, Wx,<k}] < n\zkAxf] < K. (5.15) 

We now claim [zaJ > b Otherwise, \zx^ \ < 1. By the definition of A/, It/aJ > ^ and \yxi-i\ < I for 
i > 0. Then by Zk = yk - 9Af{yk), 

\zxf = \yxi\^-29Aylj{yx,)+ \f{yxi)\‘^A‘^9^ 

> ( 1 + 76lA)|yAj^ 

> {l + -f9A)f>f, 

which is a contradiction. By (I5.15p . for any k > 0 

PIA; < A:} < - 75 - —>■ 0, Z —>■ 00 , (5.16) 

r 

that is, as I —>• 00 , A; t 00 a.s. Define /ta = 7 — (1 — 29)KA — 0.5(tA. For A < A|, ^a > 0. Using 
the linear growth condition (j5.5p gives 

\zk? = \yk - 6'/(yfc)Ap < (1 + 9AxfKf\yk\‘^. 
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Substituting this into (I5.13P yields that for A < AJ 


kfc+iP < (l - 


/taA 


which implies 


(1 + baVk)‘ 


k-l 


+ mfc, 


< |zoP + 


j=0 


where 7 a = -^ log (l “ (i+gA^)^ )' R-opl^cing k hy k A Xi yields 


k-l 


g7A(fcAA0A|^^^^j2 < |^o|2 + ^e'^'^^^l[o7,](j)mj• 


Taking expectation gives 


jE[g7A(A:AA7A|^^^^j2] 


By Fatou’s Lemma, letting / ^ oo, we have 

e'y^’^^E[\zk\^] < E|zoP, 


(5.17) 


which together with (15.81) implies the desired assertion. 

Case (ii) 9 G (1/2,1]; Since 9 G (1/2,1], by Af{yk) = (yfc — '^fc)/^, we obtain from the definition 
of nik 

rrik = 2(^yk + ^-Y^{yk- Zk),giyk)^Wk + ^L^giyk)i\^'Wk\'^ - A)'^ 

+\9{yk)\^{\^Wk\‘^ - A) + {g{yk)Awk, L^g{yk){\Awk\‘^ - A)) 
+\\L^g{yk)mAwk\^-Af-2A‘^] 

and from (15.3p 

\g{yk)\^ < -2ylf{yk) = -^[yl{yk - Zk)] < ^(3|2/feP + \zk?)- (5.18) 

This, together with ()5.1I) . produces that there is a positive constants C_ such that 

\mk\ < 01 + lArcfcl^ + \Awk\'^]\yk?‘ + Q_[l + \Awk\^]\zk\^■ (5.19) 


We therefore have E[mfcl[o 7 (](A;)] = 0. By (I5.10p and Af{yk) = (yk - Zk)/9, 


\Zk+l\'^ < 


02 


\Zk\ + 


1 - 29 


02 


- (7 - 0.5(tA)A lyfcp + 


40-2 


02 


-Zk yk + rrik. 


(5.20) 


Note that A < ^ implies 20 — 1 + (7 — O.5itA)A02 > 0. Since 

20 - 1 + (7 - O.5aA)A02 20-1 

H Vk < - \rn —^^- \yk\ + 


29 - 1 


20 - 1 + (7 - O.5aA)A02 


\Zk\ 
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we get from (15.201) 


\Zk+l\^ < 


{i-ef 


+ 


{29 - If 


\zk\ + 


02 02(20-1 +( 7 -O.5aA)A02)J 

(1-20)(7-O.5ctA)A 


\zk?‘ + ruk 


20 - 1 + (7 - O.5fjA)A02 


\zk\ +mk, 


(5.21) 


By the similar techniques used in case of 0 G [0,1/2], we can prove that Xi '[ 00 a.s., as I —>■ 00 . 
Note that (|5.21l) implies 


k-l 

1=0 

where 7 a = —-^log ^1 — obtain the desired assertion by repeating 

the proof process of the case (i). □ 

Remark 5.1. Theorem 15.21 shows that the two theta-Milstein schemes can share the exponential 
mean-square stability of the exact solution. In fact, the upper bound 7 of the decay rate can also 
be reproduced arbitrarily accurately for sufficiently small stepsize A since limA^.o 7 A = 7 - 

Let us now examine the linear scalar system 


dx{t) = fix{t)dt + cx{t)dw{t) 


(5.22) 


It is known that mean-square stability for (|5.22p is equivalent to 

2// + < 0. (5.23) 

Moreover, it is easy to observe that the two classes of Milstein schemes are equivalents for the linear 
SDE (I5.22p . Here we rewrite this scheme as follows: 

c 2 

Vk+I = yk + O^^Vk+i + (1 - 0)y^yk + cykAwk + —yk[\Awkf - A]. (5.24) 

For the linear scalar SDE (I5.22p . we have the following stability theorem. 

Theorem 5.3. Let condition (|5.23l) hold. If 


{2n + c^) + ^c^A + (1 - 29)^^A < 0, 
then the Milstein scheme ()5.24p holds 


where 7 a 


nykf = e-^^^^nxo\f 

1 / [2/i + c2 + (l-20)/x2A + l/2c^A]A 

+- (1 - - 


(5.25) 
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Proof. Rearranging equation (|5.24l) gives 


1 / \ 

Vk+i = “ 9)nA + cAwk + —[\Awk\^ - A]jyk. 

Note that yk is ^kA- measurable and Awk is independent of dkA- We therefore have 

[l + (l-0)/rA]2 + c2A + 4A2 


lE|yfc+ 


i| — 


1 + 


(1 - 

[2/i + c2 + (1 - 26i)p2A + l/2c^A]A 


(1-0^A)2 






a^ol 


where we used KAwk = 0, K\Awk\'^ = A, E(Ar(;fc)^ = 0 and E|Art;fc|^ = 3A^. □ 

Remark 5.2. Theorem 15.31 shows that (i) for 6 G [0,1/2], if A < A* := q 5 c 4 ]^(i~ 20 )/^^ ’ Milstein 
scheme (|5.24p shares the exponential mean-square stability of the exact solution, (ii) For 9 G 
(1/2,1], if jjL^ < 2 ( 29 - 1 ) diffusion term plays a crucial role), then for A < A*, Milstein scheme 
(j5.24p is exponentially mean-square stable, and if > 2 ( 29 - 1 ) drift term plays a crucial role), 
then Milstein scheme ()5.24|) is exponentially mean-square stable unconditionally. These results are 
coincident with Theorem 2.1 in [35]. Theorem 15.31 also presents the exponential decay rate 7 a of 
theta-Milstein scheme. 
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